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Summary  -  The  Marciniak  and  Kuczynski  (M-K)  analysis,  which  associates 
localized  necking  in  sheet  metal  with  pre-existing  imperfections,  has  been 
modified  such  that  the  imperfection  or  groove  is  inclined  (rather  than 
normal)  to  the  major  principal  strain  axes.  The  analysis,  which  is  closely 
related  to  the  calculations  of  Hutchinson  and  Neale,  takes  into  account  the 
requirement  for  localized  necking  that  a  direction  of  zero-extension  is 
inclined  for  strain  states  between  uniaxial  tension  and  plane  strain. 

Based  on  the  inclined  imperfection  model,  calculations  are  presented  for 
the  negative  minor  strain  regime  of  a  forming  limit  diagram  showing  the 
dependence  of  localized  necking  on  work  hardening,  strain  rate  sensitivity, 
plastic  anisotropy,  and  imperfection  factor.  Theoretical  calculations  of 
forming  limit  curves  indicate  good  agreement  between  the  model  and  previously 
obtained  experimental  data  for  several  alloys.  Such  agreement  is  obtained 
for  realistically  small  imperfections,  although  the  limit  strains  predicted 
by  the  model  itself  are  relatively  insensitive  to  the  magnitude  of  the 
imperfection. 
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NOTATION 


work  hardening  exponent  in  the  Swift  equation 

strain  rate  sensitivity 

plastic  anisotropy  parameter 

imperfection  factor 

constant  in  the  Swift  equation 

strength  constant  in  the  Swift  equation 

thickness  of  sheet  in  reagion  A  and  B  respectively 

major  and  minor  limit  strain 

the  angle  between  the  normal  of  the  imperfection  and  the  axis 

Hill's  angle  of  direction  of  zero-extension 

stress  ratio  =  c^/o^ 

strain  ratio  = 

effective  stress  and  strain 

principal  stresses 

principal  strains 

stress  components  in  the  x-y  coordinate  system 
strain  components  in  the  x-y  coordinate  system 


Introduction 


Sheet  materials  deforming  under  multiaxial  states  of  stress,  as  in  sheet 
metal  forming  operations,  usually  fail  by  localized  necking.  The  current 
interest  in  understanding  sheet  metal  formability  has  led  to  several  theoretical 
analyses  of  localized  necking  based  on  different  criteria.  These  localized 
necking  criteria  include:  a  localized  shear  zone  along  a  direction  of  zero- 
extension  [1],  material  izr$>erfections  12] ,  or  the  presence  of  a  vertex  on 
the  yield  surface  [3] ,  and  void  growth  [4] . 

Localized  necking  along  a  direction  of  zero-extension  was  originally 
proposed  by  Hill  [1] .  Hill's  theory  predicts  that  the  maximum  principal 
strain  ej  prior  to  localized  necking  (i.e.,  the  limit  strain)  has  a  magnitude 
of  e*  =  n  at  plane  strain  and  increases  to  £*  =  (1+R)n  for  the  uniaxial  tension 
deformation  of  sheet  exhibiting  normal  anisotropy  with  a  plastic  anisotropy 
parameter  R,  which  is  defined  as  the  ratio  of  the  width  strain  to  thickness 
strain  of  sheet  specimens  deformed  under  uniaxial  tension.  For  plastically 
isotropic  material  (R=l) ,  the  limit  strain  at  uniaxial  tension  thus  reduces 
to  the  well-known  e*  =  2n  expression.  Hill's  theory,  however,  cannot  explain 
the  phenomenon  of  localized  necking  under  biaxial  states  of  stress.  A  direction 
of  zero-extension  does  not  exist  when  both  and  e2  are  greater  than  zero 
and  thus  Hill's  theory  is  limited  to  the  negative,  minor  strain  regime  of  a 
forming  limit  diagram  (p  <  0) . 

Strain  localization  developed  by  local  weakness  of  material  was  first 
proposed  by  Marciniak  and  Kuczynski  (M-K)  [21  and  extended  by  Sowerby  and 
Duncan  [5].  In  the  M-K  analysis,  a  material  imperfection  in  the  form  of  a 
groove  or  trough  is  postulated  to  occur  perpendicular  to  the  direction  of  major 
principal  strain  (s^) .  Imposing  the  same  e2  inside  and  outside  the  groove  while 
proportional  straining  is  maintained  outside  the  groove,  M-K  have  shown  that 
deformation  within  the  groove  occurs  at  a  faster  rate  than  the  rest  of  the 
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sheet.  The  concentration  of  strain  (e^)  within  the  groove  eventually  leads 
to  the  plane  strain  condition  (de2  "  0)  within  the  groove  and  to  localized 
necking.  The  M-K  model  is  thus  able  to  explain  localized  necking  under  biaxial 
stresses.  Recently ,  Ghosh  [6]  used  the  M-K  model  to  calculate  the  forming  limit 
diagrams  for  A-K  steel  and  70-30  brass.  He  found  reasonable  agreement  between 
the  calculated  and  esq>erimental  forming  limit  diagrams  of  these  alloys.  Thus, 
he  concluded  that  Marciniak  and  Kuczynski's  model  of  imperfection  growth  can 
be  used  to  predict  sheet  metal  forming  limits,  provided  that  an  accurate 
constitutive  equation  is  used  and  imperfections  of  metallurgical  origin  are 
accounted  for. 

An  attractive  feature  of  the  M-K  analysis  is  that  the  concept  of  the  im¬ 
perfection  allows  one  to  interpret  both  the  negative  and  positive  sides  of  a 
forming  limit  diagram  (-  <  p  <  1)  in  terms  of  a  single  criteria  (as  opposed 

to  relying  on  Hill's  theory  for  -  <  p  <  0  and  another  criteria  for  biaxial 

tension).  In  the  M-K  analysis,  localized  necking  is  developed  from  the  strain 
concentration  within  a  material  imperfection  in  the  form  of  a  groove  which  is 
assumed  to  align  perpendicular  to  the  major  principal  strain  axis  (£i) .  In 
biaxial  straining  of  sheet  metal  (£i  >  ^2  >  °) /-  localized  necking  does  occur 
normal  to  the  major  principal  strain  axis,  as  assumed  in  the  M-K  analysis. 

However,  localized  necking  usually  occurs  at  an  angle  inclined  to  the  major 
principal  strain  axis  when  the  loading  places  the  strain  path  in  the  negative 
minor  strain  regime  of  the  forming  limit  diagram  (i.e.,  p  <  0).  The  angle  of 
inclination  depends  on  the  state  of  stress  as  well  as  plastic  anisotropy  and 
thus  the  R-value  of  the  material.  For  a  plastically  isotropic  material  in 
uniaxial  tension,  localized  necking  occurs  inclined  to  the  stress  axis  at  an 
angle  of  about  55°,  which  differs  considerably  from  the  90°  assumed  by  M-K. 

Localized  necking  involving  an  inclined  groove  when  p  <  0  has  previously 


been  undertaken  by  Hutchinson  and  Neale  [7].  In  their  analysis,  a  groove  is 


is  initially  inclined  at  an  angle  to  the  0i  (or  £i)  axis*  During  straining, 
the  groove  is  rotated  toward  the  axis.  Their  calculations  show  that  the 
limit  strain  is  influenced  by  the  orientation  of  the  groove  with  respect  to 
the  Oi  axis.  Specifically,  the  limit  strain  depends  on  the  initial  as  well  as 
the  final  groove  orientation,  and  the  limit  strain  has  thus  to  be  minimized 
with  respect  to  the  groove  orientation.  The  minimum  limit  strain  occurs  when 
the  final  orientation  of  the  groove  coincides  with  the  Hill  angle  [1]  provided 
that  the  von  Mises  flow  theory  of  plasticity  is  assumed  [7]. 

This  paper  suggests  an  alternative  and  efficient  method  of  predicting 
the  limit  strain  of  sheet  metal  failed  by  localized  necking  inclined  to  the 
axis  for  strain  paths  p  <  0.  The  formulation  of  the  present  analysis  is 
similar  to  that  of  Hutchinson  and  Neale  [71  for  the  case  of  J  flow  theory  of 
plasticity  and  therefore  a  smooth  yield  surface.  In  this  anlaysis  the  orien¬ 
tation  of  the  inclined  groove  is,  however,  rigidly  fixed  along  Hill*s  angle 
of  zero-extension.  This  simplifying  assumption  allows  the  calculation  of 
the  minimum  forming  limit  strain  without  minimizing  the  limit  strain  as  a 
function  of  the  groove  orientation.  This  permits  a  straightforward  examination 
of  the  dependence  of  localized  necking  (for  p  <  0)  on  such  factors  as  work 
hardening,  strain  rate  sensitivity,  plastic  anisotropy  and  the  size  of  the 
imperfection.  Comparison  of  analysis  with  experimental  data  can  also  be 
readily  determined. 

In  reality,  material  imperfections  are  not  likely  to  be  conveniently 
aligned  along  a  direction  of  zero-extension.  However,  we  suggest  that  localized 
necking  is  initiated  at  some  arbitrary  inclined  angle  and  quickly  rotates  into 
the  zero-extension  configuration  during  subsequent  straining.  This  analysis 
is  thus  an  approximate  solution,  but  it  does  give  predictions  of  limit  strains 
virtually  identical  to  those  calculated  by  Hutchinson  and  Neale  [7]  who  allow 


the  imperfection  to  rotate. 


The  present  analysis  relates  to  the  original  M-K  analysis  [2,5]  in  that 
it  converges  with  the  M-K  predictions  at  plane  strain.  For  deformation  under 


biaxial  tension  (£i  >  £2  >  0) ,  the  orientation  of  the  localized  neck  is  not 
inclined  but  normal  to  the  £i  axis.  Thus  previous  calculations  [2,5]  and 
comparison  [6]  with  the  positive  side  of  experimental  forming  limit  diagrams 
apply. 

Analysis 

The  present  approach  follows  the  procedure  used  by  Marciniak,  Kuczynski 
and  Pokora  [8] .  The  analysis  is  within  the  framework  of  generalized  plane  stress 
and  is  a  "long- wave length”  approximation  as  discussed  by  Hutchinson,  Neale  and 
Needleman  [9].*  As  shown  in  Fig.  1,  a  material  imperfection  is  postulated 
to  occur  on  an  angle  inclined  to  the  principal  stress  axis  (0^)  of  a  sheet 
metal  under  biaxial  stresses  of  0^  and  02-  The  inclined  imperfection  is 
designated  as  region  B  and  its  thickness  is  tfi.  The  region  outside  the  imper¬ 
fection  is  considered  to  be  uniform  and  is  referred  as  region  A;  its  thickness 
is  designated  as  tA.  The  angle  between  the  normal  of  the  imperfection  and  the 
0^  axis  is  noted  as  <p . 

Following  Marciniak  et  al.'s  approach,  [8]  the  equilibrium  requirement 
between  regions  A  and  B  gives : 

H  °XB  =  tPi  axA 

Eqn.  (1) 

and  fcB  °xyB  =  fcA  °xyA* 

where  0V  ,  0V  ,  0W  and  0W  are  the  respective  stress  components  in  regions 
XA  B  357  A 

A  and  B.  The  stresses  0V  ,  0„  and  0W  are  obtained  by  stress  transformation 

A  yA  *yA 

of  0^  and  02  into  the  x-y  coordinate  system. 


*The  basis  for  a  long  wavelength  approximation  is  that  the  ratio  of  the  width 
of  the  localized  neck  to  sheet  thickness  is  large  (>4)  such  that  no  stress  or 
strain  gradient  exists  across  the  thickness  of  the  sheet  [9] . 
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^  sin2  $  +  a2  cos2  <J>  Eqn.  (2) 

'A 

rxy  =  ^al"a2^  sin<J)  cos<J> 

To  relate  stresses  to  strains,  a  power  law  type  effective  stress-effective 
strain  constitutive  equation  is  used  similar  to  that  proposed  by  Swift  [10] : 

0  =  K  (e0  +  e)n  e  m  Eqn.  (3) 

where  0  =  effective  stress,  £  =  effective  strain,  £Q  -  a  prestrain  term, 
n  =  strain-hardening  exponent,  K  =  strength  constant,  £  =  effective  strain  rate 
and  m  =  strain-rate  sensitivity  exponent.  The  imperfection  is  assumed  to  be  a 
variation  in  K  x  thickness.  The  imperfection  factor,  f,  is  thus  defined  as 
(KQt3)/(K^tA) ,  where  and  KB  are  the  strength  constants  in  regions  A  and  B 
respectively.  Combining  Eqns.  (1)  and  (3),  w*  obtain: 


«V°A>  (deA/deB)m 
(0xb/°B^ 


=  f 


eo  +  CB 
+  EA 


exp[£ZB  -  eZA] 


Eqn.  (4) 


where  e_  and  e_  are  the  thickness  strains  in  regions  A  and  B  respectively. 

ZA  ZB 

The  sheet  is  assumed  to  have  normal  anisotropy  and  we  assume  that  the 
effective  stress  is  adequately  described  by  Hill's  original  yield  function  for 
anisotropic  material : * 


-  _/3  /R+l  r 

V  2/  R+2  J  ■ 


2  2R  _  .  _  2  .  2(1+2R)  _  2 

x  ‘  in ax  °y  °y  +  i+R  axy 


Eqn.  (5) 

where  R  is  the  plastic  anisotropy  parameter.  The  effective  strain  thus  becomes 


*In  order  to  study  the  influence  of  plastic  anisotropy  on  localized  necking  and 
thus  the  forming  limit  strain,  an  anisotropic  yield  function  is  required. 

Hill  [11]  has  proposed  a  new  anisotropic  yield  function  which  consists  of  a 
stress  exponent  parameter,  M.  The  value  of  the  parameter,  M,  has  to  be  determined 
by  biaxial  experiments  and  is  not  available  for  most  materials.  The  earlier 
function  is  adequate  for  R  £  1. 
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ai  -  /  1  /  ijjWRL 

/  3  (/  1+2R 

and  the  associated  flow  rules  are 
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fdexz  +  dexdey  +  dey  +  dexy  Eqn.  (6) 
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3R  de 
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Eqn.  (7) 


From  the  theory  of  Marciniak  et  al  [8] ,  the  constraint  condition  between  regions 


A  and  B  is 


deyA  =  deyB 


Eqn .  (8) 


If  the  inperfection  orientation  <)>  is  taken  along  a  direction  of  zero-extension 
such  that  plane  strain  condition  is  satisfied  as  d£y  =  d€y  =  0.  The  angle  ^ 


is  designated  as  <j>*  and  is  given  by  (1): 


4>*  -  90°  -  tan-1  /  a  (  R  * 

/  ,1+R, 

/  hr*  a  ~ 1 


Eqn.  (9) 


After  a  lengthy  mathematical  manipulation ,  the  following  expression  for  describing 
the  deformation  of  sheet  metal  exhibiting  localized  necking  inclined  to  the 
axis  is  obtained: 


d eA  [e0+eA]  exp[-c£ft]  -  f[e0+eBl  exp[-F  -  eB]deB 


Eqn.  (10) 


where 
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(1-g)  sincft*  cosft* 
cos^<{>*  +  a  sin*<j>* 


(1-p)  sin<fr*  costp* 


e2 
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Results 


As  shown  in  Fig.  2,  the  limit  strains  calculated  from  Eqn.  (1),  which  are 
solved  numerically  by  a  simple  incremental  method ,  are  almost  identical  with 
those  previously  obtained  in  a  more  complex  and  rigorous  method  by  Hutchinson 
and  Neale  [7]  for  the  case  of  an  imperfection  which  is  allowed  to  rotate.  The 
results  are  shown  in  Fig.  2  for  the  negative  portion  of  the  FLD,  recalling  that 
inclined  localized  necking  occurs  only  when  p  <  0. 

Comparisons  of  calculated  limit  strains  between  M-K  analysis  with  perpen¬ 
dicular  groove  ((f>=0)  and  the  inclined  groove  model  have  also  been  made  and  an 
example  of  their  differences  is  depicted  in  Fig.  3  for  the  case  of  A-K  steel. 

When  p  =  ^2/el  =  0  as  in  plane  strain  condition,  Eqn.  (9)  gives  4>*  =  0  and 
therefore  the  inclined  imperfection  model  is  equivalent  to  the  M-K  analysis  and 
both  models  given  the  same  effective  limit  strain  for  plane  strain  loading.  The 
M-K  model,  however,  predicts  a  much  greater  effective  limit  strain  as  the  strain 
ratio  p  becomes  more  negative  for  the  small  irrperfection  (f  =  .998)  used.  For 
example,  in  uniaxial  tension  (p  =  -0.60),  M-K  predicts  a  limit  strain  of  G|  =  2.8 
while  the  present  analysis  predicts  G*  =  0.74.  Experimentally,  the  limit  strain 
for  localized  necking  of  A-K  steel  is  G*  =  0.79  [12).  The  M-K  analysis  can  be 
brought  into  better  agreement  with  the  experimental  value  by  using  a  larger 
value  of  f  (f  <  .98),  but  this  seems  unrealistically  large.  In  addition, 
localized  necking  in  A-K  steel  is  observed  to  be  inclined  at  an  angle  of  <J>  =  27° 
which  compares  favorably  to  <J>  =  37°  from  Eqnf  (9)  . 

The  distribution  of  the  strain  components  of  the  effective  strain  attained 

outside  the  imperfection  site  at  the  onset  of  plastic  instability  is  depicted 

_  * 

in  Fig.  4.  The  limit  effective  strain  eA  ,  computed  by  the  inclined  imperfection 
model,  is  seen  to  increase  as  the  strain  state  (p)  becomes  more  negative.  This 
trend,  however,  is  not  always  observed  as  far  as  individual  strain  components 
are  concerned.  The  strain  components  £*  and  £*  are  both  constant  and  independent 
of  p  (£*  *  0.30  and  e*  »  0.0).  The  shear  strain  component  e*y,  on  the  other  hand. 


Fig.  2.  A  comparison  of  the  calculated  limit  strains  based  on  the 
present  analysis  and  that  of  Hutchinson  and  Neale  (71 . 
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The  distribution  of  the  strain  components  of  the  effective 
limit  strain  as  a  function  of  strain  ratio. 


increases  very  rapidly  as  p  becomes  more  negative.  In  fact,  the  increase  in 
is  largely  due  to  the  increase  in  e£y  as  the  stress  state  turns  toward  the 
uniaxial  tension  (Fig.  4) . 


On  the  Influence  of  Materials  Properties 

The  influence  of  material  properties  such  as  n,  m  and  R  on  localized 
necking  can  be  readily  demonstrated  using  the  inclined  imper feet ion  model. 

The  strain  ratio  has  been  chosen  at  uniaxial  tension  (p  =  -  j^)  and  a  very  small 
imperfection  is  used  (f=0.9999) .  Figure  5  depicts  the  influence  of  the  strain 
hardening  exponent,  n,  on  the  flow  localization  process,  while  maintaining  the 
other  parameters  constant  (m=0,  eo«0,  R^l)  .  The  results  indicate  that  the 
strain  within  the  imperfection  site  is  always  larger  than  that  outside  the 
imperfection ,  and  when  flow  instability  occurs,  strain  is  localized  within  the 
imperfection  while  the  region  outside  ceases  to  deform  at  the  limit  strain. 

Figure  5a  clearly  shows  that  a  higher  strain  hardening  exponent,  n,  can  delay 
the  onset  of  instability  to  higher  strain  value  and  thus  enhance  the  limit  strain. 
Interestingly,  the  limit  strains  calculated  for  each  of  the  n  values  is  just 
slightly  smaller  than  2n  which  is  the  prediction  of  Hill's  theory  for  isotropic 
material  with  R=l. 

The  influence  of  strain  rate  hardening  exponent,  m,  on  the  imperfection 
growth  (n=0.20  in  this  case)  is  shown  in  Fig. 5b.  The  strain  localization  process 
occurs  quite  abruptly  at  m=0.  However,  strain  rate  hardening  appears  to  have 
a  stabilizing  effect  on  strain  localization  as  it  is  delayed  and  occurs  at  a 
much  more  gradual  rate  as  m  increases.  The  consequence  of  this  stabilizing 
effect  is  the  enhancement  of  the  limit  strain,  which  increases  from  slightly 
less  than  2n  (ej*0.38)  at  m  *  0.0  to  much  greater  than  2n  (ej«0.84)  at  m  *  0.02. 

The  plastic  anisotropy  parameter,  R,  can  also  exert  a  beneficial  influence 
by  delaying  localized  neck  formation  for  p  <  0.  Figure 5c  shows  that  the  limit 


14 


strain  is  enhanced  significantly  as  the  R  value  increases.  The  dependence  of  the 
limit  strain  on  the  R-value  is  also  predicted  by  Hill's  theory.  The  difference 
between  the  predicted  effect  of  R  according  to  Hill's  theory  vs.  the  inclined 
imperfection  model  is  small  when  m  =  0.0  but  the  difference  becomes  greater  as 
m  increases.  Thus,  for  moderate  m  (.01)  values,  the  present  analysis  predicts 
a  limit  strain  in  uniaxial  tension  which  is  about  60%  greater  than  that  predicted 
by  Hill's  theory. 

The  onset  of  localized  necking  is  obviously  influenced  by  the  severity 
of  the  imperfection.  Figure  6  depicts  the  dependence  of  the  limit  strain  on  the 
imperfection  factor  f.*  Similar  to  the  M-K  analysis,  the  inclined  imperfection 
model  is  affected  by  the  imperfection  size.  The  limit  strain  €\  is  reduced  as 
the  imperfection  size  gets  larger  (f  becomes  smaller) ,  indicating  that  imper¬ 
fection  growth  is  accelerated  by  an  increasing  imperfection  size.  However, 
unlike  the  M-K  analysis,  the  inclined  imperfection  model  is  much  less  sensitive 
to  the  imperfection  factor  f.  At  a  constant  f,  the  M-K  analysis,  again,  predicts 
a  limit  strain  larger  than  that  obtained  by  the  inclined  imperfection  model. 

The  discrepancy  becomes  larger  when  the  imperfection  size  becomes  small,  i.e.r 
when  f  approaches  to  unity. 

Limit  Strain  at  Uniaxial  Tension 

Incorporating  material  parameters  of  n,  m  and  R  into  the  inclined  imper¬ 
fection  model,  the  uniaxial  tensile  limit  strains  (ej  at  a  *  0)  of  various  sheet 
metals  have  been  calculated.  The  materials  investigated  include  several  steels, 
aluminum  alloys,  brass,  and  commercially  pure  (CP)  titanium.  In  all  computations, 
a  small  imperfection  of  f  *  0.998  is  used  and  ec  is  set  equal  to  zero, 

which  zero  is  equivalent  to  assuming  that  the  material  obeys  a  power  law 

•  -n-m 

hardening  behavior:  O  ■  ke  e  .  The  orientation  of  the  inperfection  is  chosen 
to  lie  along  a  direction  of  zero-extension  and  is  given  by  Eqn.  (9).  The 

♦Recall  that  f  =  kgt-g/kAtA.  The  term  f  can  be  a  geometric  groove  (in  which  case 
f  -  tg /t^)  or  a  material  imperfection  (in  which  case  f  -  kg/k^) . 
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calculated  limit  strain  is  compared  with  experimental  data  from  literature 
[5,13,14,15,16]  and  Hill's  theory  which  predicts  e*  *  (1+R)n.  All  results  cure 
summarized  and  shown  in  Table  I  together  with  the  appropriate  material  parameters . 
Ihe  uniaxial  tensile  limit  strains  calculated  from  the  inclined  imperfection  in 
most  cases  are  in  good  agreement  with  the  experimental  values  (with  the  exception 
of  1/4  hard  Al) .  Hill's  theory,  on  the  other  hand,  is  only  in  fair  agreement 
with  the  experimental  data.  The  results  indicate  that  Hill's  theory  tends  to 
overestimate  the  limit  strain  when  m=0  and  underestimates  the  limit  strain  when 
m  >  0.  The  M-K  analysis  has  also  been  used  to  calculate  the  tensile  limit 
strain  for  A-K  steel,  HSLA  steel,  2036-Al,  70-30  brass  and  CP  Ti.  The  results 
indicate  that  the  M-K  analysis  greatly  overestimates  the  tensile  limit  strains 
of  these  materials  (see  Table  *V.  A  large  imperfection  (for  example,  f  ^  .98) 
usually  can  bring  the  M-K  prediction  into  better  agreement  with  experimental 
results,  but  such  a  larqe  ip%er feet ion  is  unrealistic. 

Forming  Limit  Diagram 

The  negative  side  of  the  forming  limit  diagram  has  been  calculated  for 
A-K  steel,  2036-Al,  HSLA  steel,  70-30  steel,  and  CP  Ti  and  are  shown  in  Fig.  7a, 
b,  c,  d  and  e  respectively.  The  inperfection  sizes  are  realistically  small, 
ranging  from  0.998  to  0.99,  depending  on  the  materials.  The  results  are 
compared  with  their  counterparts  using  the  M-K  model.  Hill's  theory  and 
experimental  forming  limit  curves.  As  illustrated  in  Fig.  7a,  b,  c,  d,  and  e, 
the  forming  limit  curve  calculated  from  the  inclined  imperfection  is  in  good 
agreement  with  the  experimental  forming  limit  curve  in  all  cases  (A-K  steel, 

HSLA  steel,  2036-Al,  70-30  brass  and  CP  Ti) .  The  agreement  between  Hill's  theory 
and  observed  behavior  is  also  good  and  compares  favorably  to  that  between  the 
present  analysis  and  the  data.  M-K  analysis  with  a  perpendicular  groove  shows 
poor  agreement  with  the  experimental  results  in  all  cases.  As  before,  the 
JMK  analysis  could  be  forced  into  a  much  better  fit  with  a  larger  inperfection 
factor,  which  in  most  cases  would  be  unrealistically  large. 


17 


•O I  -0.4  -03  -O.t  -OJ  • 


-t.O  -0  4  -0.4  -0.4  -O.t  0 


«; 


Fig .  7 .  Comparisons  of  the  negative 

side  of  the  FLD  between  Hill's 
theory  [1] ,  the  inclined  imper¬ 
fection  model,  the  M-K  analysis 
and  experimental  results  for: 

(a)  A-K  steel  [12] , 

(b)  HSLA  steel  [14], 

(c)  2036  Al  [5] , 

(d)  CP  Ti  [16],  and 

(e)  70-30  brass  [5] . 


Table  1.  Theoretical  predictions  of  uniaxial  tensile  limit  strain  by  the  inclined  imperfection 
model  and  comparisons  of  the  calculated  values  with  previously  reported  experimental 
results,  Hill's  theory  and  the  M-K  analysis. 
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Discussion 

In  this  analysis  localized  necking  for  -  <  P  <  0  is  assumed  to  occur  at 

Hill's  angle  of  direction  of  zero-extension.  By  rigidly  maintaining  an  imper¬ 
fection  under  such  configuration  during  straining,  a  simple  equation  (Eqn.  10) 
has  been  obtained  for  predicting  the  limit  strains  and  thus  FLDs  involving  negative 
Z2  strains.  Good  agreement  is  obtained  by  the  present  analysis  both  with  previously 
obtained  experimental  results  and  with  the  calculations  of  Hutchinson  and  Neale  [7] . 
This  suggests  that  localized  necking,  when  induced  by  material  imperfection ,  is 
quickly  controlled  by  the  behavior  of  an  imperfection  near  the  orientation  of 
zero-extension,  i.e,,  the  bifurcation  orientation. 

A  consequence  of  restricting  the  imperfection  along  a  direction  of  zero- 
extension  is  that  the  strain  along  the  imperfection ,  £y,  must  always  be  maintained 
at  zero  (plane  strain) ,  regardless  of  the  strain  ratio  (p) .  The  state  of  stress 
of  the  imperfection  can  thus  be  viewed  as  that  of  plane  strain  superimposed  with 
a  shear  stress  T  whose  magnitude  depends  on  p  (see  Fig.  1).  Specifically,  T 
increases  as  p  becomes  more  negative,  i.e.,  as  the  imperfection  is  inclined  at 
smaller  angles  to  the  axis.  The  principal  strain  normal  to  the  imperfection 
is,  however,  limited  to  a  constant  value  (£*  *  0.3  in  Fig.  4),  independent  of  p. 

On  the  other  hand,  the  shear  strain  term  e*  becomes  more  important  as  the  strain 

xy 

ratio  becomes  more  negative.  As  evidenced  in  Fig.  4,  it  is  the  dominant  shear 

strain  term  £*  which  is  responsible  for  the  increase  of  the  effective  limit 
xy 

strain  as  p  is  decreased  from  zero  to  more  negative  values . 

In  the  negative  minor  strain  region  (p  <  0) ,  the  M-K  model  predicts  an  eA 
larger  than  that  of  the  inclined  imperfection  (Fig.  3).  According  to  Marciniak  117], 
this  is  due  to  the  fact  that  a  change  of  stress  ratio  occurs  in  the  M-K  analysis 
as  the  strain  path  within  the  imperfection  tends  toward  plane  strain.  The  result 
is  an  increase  in  the  apparent  strength  of  the  imperfection  when  Hill's  anisotropic 
yield  function  is  used.  Thus,  the  change  of  stress  ratio  within  the  imperfection 
can  stabilize  the  process  of  straining,  delay  flow  localization  and  cause  an 
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increase  of  the  effective  limit  strain  e_.  When  the  orientation  of  the 

A 

imperfection  is  at  Hill's  angle  of  a  direction  of  zero-extension,  an  initial  strain 
path  of  plane  strain  is  specified  within  the  imperfection  and  no  change  in  the 
stress  ratio  or  increase  in  apparent  strength  occurs.  As  a  result,  the  inclined 
imperfection  predicts  an  effective  limit  strain  considerably  lower  than  that 
of  the  M-K  analysis.  The  difference  is  largest  at  uniaxial  tension  and  diminishes 
to  zero  at  plane  strain  when  p  =  0. 

Influence  of  Materials  Properties  on  the  Forming  Limit 

The  influence  of  work  hardening  in  rate  and  strain  rate  sensitivity  m  on 
forming  limit  strain  is  well  known.  Materials  with  high  n  and/or  m  values  tend 
to  have  superior  sheet  formability  because  strain  is  distributed  more  uniforml  y [5] 
Figure  5a  clearly  depicts  the  beneficial  effect  of  n  on  the  limit  strain  while 
Fig. 5b  indicates  that  even  a  small  m  tends  to  stabilize  the  localized  neck  and 
enhance  the  limit  strain  by  allowing  the  material  to  deform  in  a  quasi -stable 
manner.  The  stabilizing  effect  of  m  on  localized  necking  has  also  been  discussed 
by  Ghosh  [5].  An  interesting  observation  of  Figs. 5a  and  5b is  that  the  inclined 
imperfection  model  gives  a  limit  strain  close  to  that  predicted  by  Hill's  theory 
when  the  strain  rate  effect  is  not  taken  into  account  (m=0) .  Large  discrepancies, 
however,  occur  between  the  two  analyses  when  the  strain  rate  effect  is  also 
considered.  The  inclined  imperfection  predicts  a  limit  strain  considerably  higher 
than  that  of  Hill's  theory  when  m  >  0.  Close  agreement  between  the  inclined 
imperfection  model  and  Hill’s  theory  also  exist  in  predicting  that  an  increase 
in  R  will  increase  limit  strains,  as  shown  in  Fig.  5c.  it  should  be  noted  that 
Fig. 5c  assumes  m=0  for  m  >  0,  the  inclined  imperfection  model  predicts  larger 
limit  strain  than  does  Hill's  theory. 

The  result  shown  in  Fig.  6  is  indicative  of  of  the  detrimental  effect  of 
the  imperfection  size  on  the  forming  limit  strain.  Compared  to  the  M-K  model, 
the  inclined  imperfection  is  found  to  be  much  less  sensitive  to  the  imperfection 
size.  This  is  probably  due  to  the  fact  that  in  the  M-K  analysis  a  large  imper- 
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fection  size  is  more  effective  in  changing  the  strain  path  toward  plane  strain. 
Since  an  inclined  groove  is  always  maintained  at  plane  strain,  it  is  therefore 
less  sensitive  to  the  imperfection  size,  as  illustrated  in  Fig.  8. 

Uniaxial  Tensile  Behavior  and  Forming  Limit  Diagrams 

The  limit  strain  of  various  sheet  metals  under  uniaxial  tension  loading 
can  be  reasonably  predicted  by  the  inclined  imperfection  method  when  proper 
material  parameters  of  n,  m,  and  R  are  accounted  for  an  a  realistic  characteristic 
imperfection  level  (f=0.998)  is  used  (see  Table  I).  The  good  agreement  between 
theoretical  calculations  and  experimental  results  suggests  that  localized 
necking  can  indeed  originate  from  material  imperfection.  Strictly  speaking, 
according  to  Hill’s  theory,  an  imperfection  is  not  required  for  localized  necking 
to  occur  on  sheet  metal  under  uniaxial  tension.  However,  Hill's  theory  is  not 
able  to  predict  accurately  the  limit  strain  in  cases  where  the  flow  stabilizing 
effect  of  strain  rate  hardening  is  important  (see  Fig.  5b)  .  The  M-K  model  with 
a  perpendicular  groove,  on  the  other  hand,  predicts  an  unrealistically  large 
limit  strain  at  uniaxial  tension  as  shown  in  Table  I.  By  restricting  the  groove 
to  align  along  a  direction  of  zero-extension,  the  inclined  imperfection  can 
predict  a  reasonable  uniaxial  tensile  limit  strain  which  compares  favorably 
with  the  experimental  one.  The  analysis  is  also  consistent  with  the  fact  that 
most  sheet  metals  under  uniaxial  tensile  loading  show  a  localized  neck  inclined 

j 

to  the  stress  axis. 

In  addition  to  the  fairly  good  prediction  of  localized  necking  in 
uniaxial  tension,  the  inclined  imperfection  model  predicts  with  reasonable  accuracy 
the  negative  minor  strain  regime  of  the  forming  limit  diagram  for  A-K  steel, 

HSLA  steel,  2031-A1,  70-30  brass  and  CP  Ti.  These  materials  exhibit  a  fairly 
wide  range  of  material  parameters.  The  good  agreement  between  the  present 
analysis  and  the  experimental  result  suggests  that  flow  localization  within 
material  imperfection  can  indeed  result  in  localized  necking.  Hill's  theory 


also  provides  a  fair  prediction  of  the  forming  limit  strains,  but  it  tends  to 
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give  an  underestimation  when  m  >  0  but  an  overestimation  when  m  =  0;  this  is 
principally  due  to  the  beneficial  effect  of  strain  rate  hardening  in  delaying 
flow  localization.  The  M-K  analysis,  on  the  other  hand,  takes  a  strain  path 
different  than  that  of  the  inclined  imperfection  model  by  using  a  perpendicular 
groove  and  predicts  an  unrealistically  large  forming  limit  strain. 

With  regard  to  the  entire  forming  diagram  including  the  biaxial  side,  the 
use  of  a  single  criterion,  a  small  imperfection,  is  appealing.  However,  while 
analyses  based  on  an  inclined  imperfection  (f  ^  .998)  yield  good  predictions 
of  localized  necking  for  p  <  0,  the  same  small  imperfection  factor  predicts 
limit  strains  which  are  much  greater  than  those  experimentally  observed  in 
biaxial  tension  when  p  =  1  [6] .  Possible  reasons  for  the  discrepancies  between 
experiment  and  the  M-K  predictions  at  p  =  1  have  been  discussed  elsewhere  [6] , 
and  include  the  intervention  of  fracture  and  strain-associated  changes  in  biaxial 
work  hardening  behavior.  Analysis  of  behavior  under  balanced  biaxial  tension 
thus  remains  a  problem. 

Conclusion 

The  Marciniak-Kuczynski  (M-K)  analysis  of  localized  necking,  which  is 
based  on  an  imperfection  normal  to  the  major  principal  strain  axis,  has  been 
modified  for  p  <  0  to  take  into  account  a  groove  inclined  to  the  ei  axis.  The 
orientation  of  the  inclined  imperfection  is  assumed  to  lie  along  a  direction  of 
zero-extension.  Such  an  imperfection  orientation  coincides  with  the  minimum  value 
of  the  limit  strain  obtained  by  Hutchinson  and  Neale  by  minimizing  the  limit 
strain  with  respect  to  the  orientation  of  the  imper feet ion.  Calculations  of 
forming  limits  predicted  by  the  inclined  groove  model  indicate  that,  for  a  given 
strain  state,  the  limit  strain  increases  with  increasing  work  hardening  rate  n, 
the  strain  rate  hardening  exponent  m,  the  plastic  anisotropy  parameter  R,  but 
decreased  with  increasing  imperfection  size  (i.e.,  as  f  becomes  waller).  The 
inclined  imperfections  model  has  been  found  to  be  much  less  sensitive  to  the 
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imperfection  factor,  f,  than  that  of  the  MHC  analysis,  When  applied  to  the 
calculation  of  the  uniaxial  and  forming  limit  strains  of  various  steels , 
aluminum  alloys,  70-30  brass  and  CP  Ti  in  the  negative  minor  strain  regime , 
the  inclined  imperfection  model  gives  better  agreement  with  the  experimental 
results  than  Hill's  theory  or  the  M-K  analysis.  Especially  for  the  case  of 
uniaxial  tension.  Hill's  theory  tends  to  overestimate  the  limit  strain  when 
m  =  0  but  underestimates  it  when  m  >  0.  Because  it  involves  a  change  in 
stress  ratio  within  the  imperfection  which  strengthen  it  during  straining 
even  for  p  <  0,  the  M~K  analysis,  on  the  other  hand,  predicts  a  limit  strain 
which  is  either  too  large  or  requires  an  imperfection  factor  that  is 
unrealistically  large. 
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